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1 Introduction 

Shear viscosity 77 is one of transport coefficients in fluid dynamics, which char- 
acterizes how viscous the system is in the presence of flow gradient [lj. Since, 
in general, shear viscosity is inversely proportional to scattering cross-section, 
strongly interacting systems have smaller viscosity than weakly interacting 
ones. Recently a universal minimum bound for the ratio of 77 to entropy den- 
sity s has been proposed by Kovtun, Son, and Starinets (KSS) [2]. The bound, 
V/ S — 2~) i s conjectured to be satisfied for a large class of strongly interact- 
ing quantum field theories whose dual descriptions in string theory involve 
black holes in anti-de Sitter space [3[ [H [5] [6] . Note that 77/s is more physical 
quantity than 77 itself because the ratio appears as a diffusion constant of fluid 
equations. 

In experiments, rj/s close to the minimum bound were found in relativis- 
tic heavy ion collisions (RHIC) [7J [H [5]. This discovery came as a surprise. 
Traditionally, quark gluon plasma (QGP), a phase of the quantum chromody- 
namics (QCD) above the deconfinement temperature T c ~ 170 MeV at zero 
baryon density [lOj , had been thought to be in weak interaction regime. Partly 
because lattice QCD simulations of the QGP equation of state above 2T C are 
not inconsistent with that of an ideal gas of massless particles, e = 3p, where 
e is the the energy density and p is the pressure of the system [10 . However, 
recent analyses of the elliptic flow generated by non-central collisions in RHIC 
[H [9] and lattice simulations of a gluon plasma [11] yielded 77/s close to the 
the minimum bound at just above T c . This suggests QGP is strongly inter- 
acting at this temperature. Other implications for strong coupling can be 
seen in sharp peaks of mesonic correlators [T2j [T3l [14] , and in discussions of 
the possible microscopic structure of such a state [16l [TjJ [TH QjH [20] . 

Given this situation, one naturally gets interested in how 77/s behaves 
as temperature approaches T c from below, supposing that 77/s of QCD was 



1 See Ref. [15] for other possibility to yield a small viscosity in expanding QGP. 
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already close to the minimum bound at just above T c , and how it relates to 
change of the effective degrees of freedom through a phase transition or cross 
over. 

To explore these issues, we use chiral perturbation theory (ChPT) and the 
linearized Boltzmann equation to perform a model independent calculation to 
the rj/s of QCD in the confinement phase. Earlier attempts to compute meson 
matter viscosity using the Boltzmann equation and phenomenological phase 
shifts in the context of RHIC hydrodynamical evolution after freeze out can 
be found in Refs. [551 H3 [55] • In the deconfincment phase, state of the art 
perturbative QCD calculations of r] can be found in Rcfs. [29l [30] , 

1.1 Shear viscosity in Fluid dynamics 

For later understanding it might be worth mentioning basic properties of shear 
viscosity. Fluid dynamics describes non-equilibrium system where evolution in 
space and time occurs in macroscopic scale, rj appears as a phenomenological 
parameter in this scale. Usually we can select suitable theory corresponding 
to the characteristic evolution scale of our interest, as listed in table [TJ 



Table 1. Hierarchy of theories in space-time scales 



Scale 


Theory 


Micro (range of interaction) 
Meso (mean-free path) 
Macro (sound wave) 


Quantum theory, e.£ 
Kinetic theory, e.f. 
Fluid dynamics, e.f 


Linear response theory 
Boltzmann eq. 
Navier-Stokes eq. with 77 



Fluid dynamics in relativistic framework is consist of two basic equations, 
that is, conservation laws of the energy- momentum T^{x) and the number 
density (or charge) n(x): 

d,T^(x) = 0, (1) 
d^n{x)V»{x) = 0, (2) 

where V^(x) is a normalized four velocity of elementary volume, V fJ, (x)V ll (x) = 
1. In what follows, we will not care the second equation in treating pion gas, 
because pion does not carry any conserved charge. The energy-momentum 
tensor in viscous system is decomposed as 

T ia ,(x)=Tf$(x)+5T ij (x) t (3) 

where the first term describes the perfect fluid (local equilibrium; frictionless 
processes) with local pressure V(x) and energy density e(x) 



T$(x) = {V(x) + e(x)} V^{x)V v {x) - V{x)5^ 



(4) 
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, and the second term corresponds to a small deviation from local equilib- 
rium, which defines shear and bulk viscosities, r\ and £, 

6T tJ {x) = -n (viVjix) + VjViix) - • V(x)\ + C%V ■ V{x). (5) 

One can see that rj (£) is in traceless (trace) part of spatial components. These 
transport coefficients, which also includes thermal conductivity in the presence 
of conserved charge, are determined by experiments, or can be derived from 
more microscopic theories in principle. In general, it becomes more difficult 
to derive them as one begins with smaller scale. 

We shall estimate ij by considering two dimensional case depicted in Fig. [I] 
where x-component of flow velocity V x (y) varies along y direction. In other 



-it 



Fig. 1. Shear flow on x-y plane. 



words, shear viscosity is an anisotropic pressure defined by the momentum 
transfer of x-component per unit time across unit area normal to y direction, 
in which the force is facing negative x direction. Since the momentum transfer 
occurs within the mean-free path A, 

rrrav dV x (y) dp x p y dV x (y) 

ST =- 11 ^y-^- nV ^ X = ~V^y-> (6) 

where we have used A = l/(na) with a and n being scattering cross-section 
and density. Thus, we can roughly estimate it by dimensional analysis up to 
a numerical factor, 

rp rp3 

where we have reduced momentum to temperature p ~ T, and T is the 
scattering amplitude. It shows that strong interaction (small A) causes quick 
and local equilibrium, implying that the system can be well described by 
hydrodynamics. In the case of weak interaction (large A), on the other hand, 
a large number of particles (thermal excitations) included in A join to relax 
the gradient of flow velocity, so that, it makes 77 finite. 



2 In general, in broken phases such as superfluid or chiral condensed phase, there 
might be an another term attributed to the superfluid mode independent from 

n(0) 



the above T^J in gj). For details, see Refs. [3Tjl32l l33 
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1.2 Shear viscosity in Kubo formula 

If one would like to evaluate rj in field theoretical way with Lagrangian, the 
linear response theory gives a formulation. Since 77 is a response of the system 
to the external stimulation 5Tij, it is given by Kubo formula 



where T lJ the spatial off-diagonal part of energy- momentum tensor, and is 
operator derived from Lagrangian of the system. This formula gives the strict 
definition of r\ derived from microscopic Lagrangian. One might think a per- 
turbative calculation of the above correlator will give an answer for r\. But this 
is not true even for weak coupling theory. Indeed, the Kubo formula involves, 
for instance in \<f> theory, an infinite number of diagrams at the leading order 
(LO) .24]. This is due to an infrared singularity (remember that rj appears in 
macroscopic scale). However, it is proven that the summation of LO diagrams 
is equivalent to solving the linearized Boltzmann equation (for mesoscopic 
scale) with temperature dependent particle masses and scattering amplitudes 
[24] , The result shows rj ~ T 3 /A 2 . This proof extended the applicable range 
of the Boltzmann equation to higher temperature, but is restricted to weak 
coupling theories. In the case we are interested (QCD with T < 140 MeV), 
the pion mean free path is always greater than the range of interaction (~ 1 
fm) by a factor of 10 3 . Thus, even though the coupling in ChPT is too strong 
to use the result of \<fi 4 theory in [24] , the temperature is still low enough that 
the use of the Boltzmann equation is justified. In the following sections, we 
will show shear viscosity of pion gas evaluated from the linearized Boltzmann 
equation to LO in ChPT, and develop arguments on results, which is all based 
on our recent paper Ref. [25] . 

2 Linearized Boltzmann Equation for Low Energy QCD 

In the hadronic phase of QCD with zero baryon-number density, the dominant 
degrees of freedom are the lightest hadrons — the pions. The pion mass m, = 
139 MeV is much lighter than the mass of the next lightest hadron — the kaon 
whose mass is 495 MeV. Given that T c is only ~ 170 MeV, it is sufficient to 
just consider the pions in the calculation of thermodynamical quantities and 
transport coefficients for T<T C . 

The interaction between pions can be described by chiral perturbation the- 
ory (ChPT) in a systematic expansion in energy and quark (u and d quark) 
masses [21j [22l [23] . ChPT is a low energy effective field theory of QCD. It de- 
scribes pions as Nambu-Goldstone bosons of the spontaneously broken chiral 
symmetry. At T <C T Cl the temperature dependence in tttt scattering can be 




(8) 
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calculated systematically. At T — T c , however, the theory breaks down due 
to the restoration of chiral symmetry|f| 

Since unlike the Kubo formula the Boltzmann equation requires semi- 
classical descriptions of particles with definite position, energy and momen- 
tum except during brief collisions, the mean-free path is required to be 
much greater than the range of interaction. Thus the Boltzmann equation 
is usually limited to low temperature systems. The Boltzmann equation 
describes the evolution of the isospin averaged pion distribution function 
/ = /(x, p,i) = f p (x) (a function of space, time and momentum) as 

CWpl?) = T / dri ^P { /l/2 ^ + W + M + + /a)/s/ P } 

A J 123 

(9) 

where E p = y^p 2 + m 2 and = 3 is the degeneracy factor for three pions, 

1 JL rf3l,. 

dr l2;3p = — \t? n j^e-) x m 

and where T is the scattering amplitude for particles with momenta 1,2 — > 
3, p. LHS of (|9|) gives temporal change of f P (x), and it should be equal to 
change by collision (RHS) in which the first (second) term corresponds to 
gain (loss) of particles of momentum p by the collision. 

In ChPT, the LO isospin averaged tttt scattering amplitude in terms of 
Mandelstam variables (s, t, and it) is given by 

|T| 2 = i £ (21 +i)|T<'>| 2 = -lj {21m* +9 S 2 -24mls + 3(t-uf} . 

^ 7=0,1,2 JJ * 

(11) 

The temperature dependence in pion mass and pion scattering amplitudes can 
be treated as higher order corrections. 

The distribution function can be decomposed into the local equilibrium 
f p °\x) and a deviation from it, 

f p (x) = f^(x) + Sf p (x), (12) 

where fp°\x) — {epi^W^i^P 11 — l) with (3(x) being the inverse temperature 
and V^(x) the four velocity at the space-time point x. A small deviation 
Sf p (x) from the local equilibrium is parameterized as 

f p (x) = /f (z) [l - {l + fW(x)} x P (x)} , (13) 
while the energy momentum tensor is 



3 The QCD chiral restoration happens almost as soon as the deconfinement tran- 
sition at zero baryon density. We do not distinguish their critical temperatures. 
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T,u(x) = T$(x) + ST^(x) =g v j [if^ + Sf p (x)] . (14) 

We will choose the ~V(x) = frame for the point x. This implies d„V° = 
after taking a derivative on V^^V^lx) = 1. Furthermore, the conservation 
law at equilibrium d^T^"]^^ = allows us to replace d t fi{x) and <9 t V(x) 
by terms proportional to V • V(x) and W(3(x). Thus, to the first order in a 
derivative expansion, Xp( x ) can be parameterized as 

— i-y = A(p)V ■ V(ai) + B ld (p) I <% 1 , (15) 

where = B{p) {pipj — and i and j are spatial indexes. 

Substituting (| 1 5|) into the Boltzmann equation, one obtains a linearized 
equation for £>(p), 

(piPj - ^ijP 2 ^ = ^ dri2 ;3 p(l + ni)(l + n 2 )n 3 (l + n p ) _1 

x [B tf (p) + Sy(fc 3 ) - By(fe) - By(fci)] 
= ^4 [B] , (16) 

where we have dropped the factor (VjVj(x) + VjVi(a;) — trace) contracting 
both sides of the equation and write ff^ (x) at this point as m — (e l3Ei — l) 1 . 
There is another integral equation associated with V- ~V(x) which is related to 
the bulk viscosity ( that will not be discussed here. The V • (3 and 9jV terms 
in p^d^fp ^ will cancel each other by the energy momentum conservation in 
local equilibrium mentioned above. 

After putting everything together and comparing the energy-momentum 
tensor |5j) in fluid dynamics and (|14j) in kinetic theory, we obtain 

U " ' ^ P 1 -n p (l + n p )B ij (p) (p^ - ^ij-p 2 ) (17) 



10 J (2n)3E p PK ■ p> '"*'Vi\W3 3 ".jj 
9lP f d 3 p 1 



10 J (2tt) 3 E, 



-n p (1 + n p ) B l3 {p)F lJ [B] ee gi{B\F[B]) . (18) 



Here one can see immediately that if the scattering cross section is scaled by 
a factor a, 

dr 12 -3p -> a {dr 12 - 3p ) , (19) 
then Eqs. ([To]) and (fT8")l imply the following scaling 



4 A non-derivative term is not allowed since f p should be reduced to fp ^ when /3 
and become independent of x. There is no term with single spatial derivative 
on f3(x) either. The only possible term (V ■ V) /3(x) vanishes in the V(x) = 
frame. 
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B l] (p)^a- 1 B lJ (p) , (20) 
ry^crV (21) 

with i] proportional to the inverse of scattering cross-section. This non- 
perturbative result is a general feature for the linearized Boltzmann equation 
with two-body elastic scattering. 

To find a solution for B(p), one can just solve Eq. (fl6)) . But here we 
follow the approach outlined in Refs. [271 [28] to assume that B(p) is a smooth 
function which can be expanded using a specific set of orthogonal polynomials 

oo 

B(p) = |p|»J>i? {r) (*(p)), (22) 

where B^ r \z) is a polynomial up to z r and b r is its coefficient. The overall 
factor |p| y will be chosen by trial and error to get the fastest convergence. 
The orthogonality condition 

(1 + n p ) \p\yB^(z)B^(z) ex 5 r , s (23) 

can be used to construct the B^ r \z) polynomials up to normalization con- 
stants. For simplicity, we will choose B^(z) = 1. 

With this expansion, the consistency condition for B(p) in Eq. (jT8j) yields 

V = g^b r L^=gl^2b r (B^\F[B^])b s , (24) 

r r,s 

where 



2 



L(T) = TeJ (^^(l + ^IPl^Wtx S , r . (25) 

Since b r is a function of m^, /„■ and T, the & r 's in Eq. (|24p are in general 
independent functions, such that = g^Y, s ( Bir) \F[ B{s) ]) b s [ 

one can show 

that this solution of b s gives a unique solution of 77] , or equivalcntly 

V£ (0) = ^E< s(r) i # [ 5(s) ]> 6 * • ( 26 ) 

s 

This will allow us to solve for the b s and obtain the shear viscosity 

ri = g„boL<®. (27) 

In the next section, we will show that this expansion converges well, such 
that one does not need to keep many terms on the right hand side of Eq. 
If only the s = term was kept, then 
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This resultant formula clearly shows that 77 does not depend on the degeneracy 
factor g n . 

The calculation of the entropy density s is more straightforward since s, 
unlike 77, does not diverge in a free theory. In ChPT, the interaction contri- 
butions are all higher order in our LO calculation. Thus we just compute the 
s for a free pion gas: 



, 2 d \ogZ 



•=~9*P W -f-, (29) 



where the partition function Z for free pions is 
logZ 1 r ^ 3 * 





up to temperature independent terms 



e-™t , (30) 



3 Results and Discussion 
3.1 Shear viscosity of Pion gas 

We present the results for rj and r]/s of QCD at zero baryon number den- 
sity. Fig. [2Ja) shows the LO ChPT result of 77 using the linearized Boltzmann 
equation. For comparison, we also added results for a it-it scattering of em- 
pirical phase shifts (PS), in which values of parameters were employed from 
Refs. 123 I2H1 131] ■ Characteristically, the PS scattering amplitude does 
not reflect the phase transition so much to be less dependent on momentum 
than ChPT at high energy region, and incorporates resonance effect which 
is absent in ChPT. The lines with circles, squares and triangles correspond 
to keeping the first one, two and three polynomials on the right hand side 
of Eq. (f2"u| , respectively. We have used y = and z(p) = |p| for ChPT and 
y = 2 and z(p) = |p| 2 for PS to construct the polynomials. The figure shows 
the polynomial expansion converges rapidly. 

As a test of the calculation, we also show the shear viscosity result for a 

constant scattering amplitude |T| 2 = ^jjr- in Fig. |2](b) (ChPT SH), which 

corresponds to the low energy limit of ChPT and can be mapped onto X(j> 
theory of Ref. [24 where the scattering amplitude is a constant |T| 2 = A 2 . In 
Xcj) 4 theory, 77 is monotonically increasing in T. If T ^> m^, 77 oc T 3 /A 2 with 
T 3 given by dimensional analysis and A~ 2 by the scaling of coupling, as shown 
in Eq. ([7]) and in Eqs. (| 19|) and ([2Tj) . In ChPT, however, r\ is not monotonic 
in T. At T m w , the scattering amplitude is close to a constant, thus ChPT 
behaves like a Xcj) 4 theory. But as temperature gets higher to T > m,, T 
oc T 2 // 2 and thus 77 oc f^/T. At what temperature the transition from r\ 
oc T 3 to r\ oc 1/T takes place depends on the detail of dynamics. In ChPT, 
this temperature is around 20 MeV. 




20 40 60 80 100 120 140 20 40 60 80 100 



T(MeV) T(MeV) 

Fig. 2. Shear viscosity as a function of temperature, (a) for the LO ChPT (ChPT) 
and for an empirical phase-shift scattering amplitude (PS) 34 . = 139 MeV and 
.U = 93 MeV. (b) for the LO ChPT in the massless limit (ChPT ML) where f v = 87 

MeV, and for a constant scattering amplitude in low-energy limit, |T| 2 = -^yf- 
(ChPT SH). l st - 3 r d show results up to the three polynomials in Eq. (1261) . 



It is also interesting to observe the result of ChPT in massless limit, 
Fig. [2Jb) (ChPT ML). 77 blows up at T = 0. It reflects the derivative cou- 
pling typical for NG modes, see Ref. [38] for rj of phonon in CFL phase. 

3.2 -q/s and error estimation 

The radius of convergence in momentum for ChPT is typically 4^/^ ~ 1 
GeV. To translate this radius into temperature, we compute the averaged 

center of mass momentum (|p|) = J (B\p 2 \F[B]) / (B\F[B]} . We found that 
for T = 120 and 140 MeV, (|p|) ~ 460 and 530 MeV < 4tt/ ot . However, ChPT 
is supposed to break down at the chiral restoration temperature (~ 170 MeV). 
Thus our LO ChPT result can only be trusted up to T ~ 120 MeV. At the 
next-to-leading order (NLO), it is known that the isoscalar tttt scattering 
length will be increased by ~ 40% [23] . This will increase the cross section by 
~ 100% and reduce 77 by ~ 50% near threshold. This is an unusually large 
NLO correction since a typical NLO correction at threshold is < 20%. The 
large chiral corrections does not persist at the higher order. At the next-to- 
next-to-leading order (NNLO), the correction is ~ 10% [23] ■ Thus, to compute 
77 to 10% accuracy, a NLO ChPT calculation is needed. 

The LO ChPT result for 77/s is shown in Fig. [3J (line with rectangles). 
The error is estimated to be ~ 50% up to 120 MeV. 77/s is monotonically 
decreasing and reaches 0.6 at T = 120 MeV. This is similar to the behavior 
in the = case (shown as the line with rectangles) where 77/s cx f 4 /T A 
with s oc T 3 from dimensional analysis and /„• = 87 MeV in the chiral limit 
[22l [23] . For comparison, we also show the result using phenomenological tt- 
tt phase shifts [34| for 77 but free pions for s. (Our result for 77 is in good 
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-CO 6 
ST 

4 
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80 100 120 
T(MeV) 

Fig. 3. Shear viscosity to entropy density ratios as functions of temperature. Line 
with circles (rectangles) is the LO ChPT result with = 139 (0) MeV and f„ = 
93(87) MeV. Line with triangles is the result using tttt phase shifts (PS). Dashed 
line is the conjectured KSS bound 1/47T ~ 0.08. 

agreement with that of [55] for T between 60 and 120 MeV. For an earlier 
calculation using the Chapman-Enskog approximation, see Ref. [35].) This 
amounts to taking into account part of the NLO irir scattering effects but 
ignoring its temperature dependence and the interaction in s. Since not all 
the NLO effects are accounted for, this rj/s is not necessarily more accurate 
than that using the LO ChPT. The comparison, however, gives us some feeling 
of the size of error for the LO result we present here. Thus, an error of ~ 100% 
at T — 120 MeV to the LO result might be more realistic. 

3.3 Relation between T c and r//s 

Naive extrapolations of the three rj/s curves show that the 1/4tt = 0.08 mini- 
mum bound conjectured from string theory might never be reached as in phase 
shift result |36[ 137] (the first scenario) , or more interestingly, be reached at 
T ~ 210 MeV, as in the LO ChPT result (the second scenario). In both scenar- 
ios, we see no sign of violation of the universal minimum bound for rj/s below 
T c . But to really make sure the bound is valid from 120 MeV to T c , a lattice 
computation as was performed to gluon plasma above T c [TT] is needed. In the 
second scenario, assuming the bound is valid for QCD, then either a phase 
transition or cross over should occur before the minimum bound is reached at 
T ~ 200 MeV. Also, in this scenario, it seems natural for rj/s to stay close to 
the minimum bound around T c as was recently found in heavy ion collisions. 

In the second scenario, one might argue that the existence of phase tran- 
sition is already known, otherwise we will not have spontaneous symmetry 
breaking and the corresponding Nambu-Goldstone boson theory at low tem- 
perature in the first place. Indeed, it is true in the case of QCD. For a sponta- 
neous symmetry breaking theory, the general feature of rj/s we see here seems 
generic. 



b\ 1 


1 i 1 i 




ps 




ChPT [M,=0 MeV] 




ChPT [M^=139 MeV] - 
1 




4 K 
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Fig. 4. A schematic of rj/s behavior around the critical temperature. 



At asymptotically high T > T c in the deconfinement phase, the weak 
coupling nature of QCD makes rj/s get higher in comparison with that near 
T c . On the other hand, at very low T <C T c in the symmetry breaking phase, 
the Nambu-Goldstone bosons are weakly interaction at low temperature, thus 
rj/s gets lower as T approaches T c from below. A phase transition should occur 
before the extrapolated rj/s curve coming from high T reaches the bound at 
T\. Similarly, a phase transition should occur before the extrapolated rj/s 
curve coming from low T reaches the bound at T2- Thus the range of phase 
transition is T\ < T c < T2, see Fig. |4] There might be a discontinuity at T c 
due to the change of effective degrees of freedom. Some relations between T c 
and the rj/s bound are also discussed in Refs. [40]. 141 j . 



3.4 Large Nf and N c limit 

It is interesting to note that the degeneracy factor g w drops out of 77 while 
the entropy s is proportional to g„, as shown in Eqs. (|26H28j) and (|29|) . respec- 
tively. This suggests the rj/s bound might be violated if a system has a large 
particle degeneracy factor [2]. For QCD, large g„ can be obtained by having 
a large number of quark flavors Nf with g n ~ Nj. However, the existence of 
confinement demands that the number of colors N c should be of order Nf to 
have a negative QCD beta function. After using / w oc \/N~ C: the combined N c 
and Nf scaling of 77/s is 

S K g-fi K Nj ' (31) 

which is of order one. Thus QCD with large N c and Nf can still be consistent 
with the rj/s bound below T c . 
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